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ON THE FINE STRUCTURE OF G-BAND BRIGHT FEATURES 
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ABSTRACT 

We present new results from first observations of the quiet solar photosphere performed through the Interfer- 
ometric BIdimensional Spectrometer (IBIS) in spectropolarimetric mode. IBIS allowed us to measure the four 
Stokes parameters in the Fe I 630. 15 nm and Fe I 630.25 nm lines with high spatial and spectral resolutions for 
53 minutes; the polarimetric sensitivity achieved by the instrument is 3 x 10~ 3 the continuum intensity level. 

We focus on the correlation which emerges between G-band bright feature brightness and magnetic filling 
factor of ~ 10 3 G (kG) fields derived by inverting Stokes / and V profiles. More in detail, we present the 
correlation first in a pixel-by-pixel study of a ~ 3 arcsec wide bright feature (a small network patch) and then 
we show that such a result can be extended to all the bright features found in the dataset at any instant of the 
time sequence. The higher the kG filling factor associated to a feature the higher the brightness of the feature 
itself. Filling factors up to ~ 35% are obtained for the brightest features. Considering the values of the filling 
factors derived from the inversion analysis of spectropolarimetric data and the brightness variation observed in 
G-band data we put forward an upper limit for the smallest scale over which magnetic flux concentrations in 
intergranular lanes produce a G-band brightness enhancement (~ 0.1")- Moreover, the brightness saturation 
observed for feature sizes comparable to the resolution of the observations is compatible with large G-band 
bright features being clusters of sub-arcsecond bright points. This conclusion deserves to be confirmed by 
forthcoming spectropolarimetric observations at higher spatial resolution. 
Subject headings: Sun: magnetic fields — Sun: photosphere — Techniques: polarimetric 



1. INTRODUCTION 

Solar magnetic fields manifest in the photosphere in a 
large variety of structures. At sub-arcsecond spatial scales 
they can appear concentrated in bright and roundish fea- 
tures, located in convective downflow regions, preferen- 
tially at the verte xes of granules . Thes e brig ht features, 
first o bserved by [Dunn & Z irkeri d 19731) and iMehltretteJ 
d 19741) in the network, are ubiquitous on the solar pho- 
tosphere dSanchez Almeida et all 12004 Ide Wiin et alj 1 20051: 
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IBovelet & Wiehri I2008t [Sanchez Alm eida et all 120 1 ob and 
their aggregation in ac tive regions form, on larger spa- 
tial scales, faculae ( e.g., iFrazier & Stenfloll 19721 iMehltretterl 
119741 iRabinl 119921) . Their photometric properties are 
largely depende nt on the spectral range in which they are 
observed (e.g.. ISutterlin et alj Il999t iLanghans et alj 12004 
iTritschler & Uitenbroekl l2006t iBeck et alj 12007b ! In par- 
ticular, they show enhanced contrast with respect to the 
quiet Sun average intensity (typically 30%) when ob- 
served in the G-band (around 430.8 n m) and are there- 
fore r eferred to as G-band bright points ( Mulle r~& Roudierl 
1984). Due to their small size, the investigation of 
these features has been long hampered by spatial resolu- 
tion limits imposed by instrumentation and atmospheric see- 
ing. During the last decade, the development of Adap- 
tive Op tic Systems as well as post-facto r estoration tech- 
nique s (lLofdahletalJ [19981: lLofdahll l200l Ivan Noort et aTJ 
120051) has allowed stu d ies that focused on the photomet- 
ric (e.g., iKellerl [1991 IBerger et all [l99l IBerger & Titlel 
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20071: Ide Wiin et alj|2008h IViticchie et alj|2009l) properties of 
these small-size magnetic flux concentrations. In particular, 
observations have r evealed that these are highly structured. 
Berge reFal] d200 4) reported that in active regions, at loca- 
tions of high magnetic flux density correspond in both Ca and 
G-band images a large variety of features that they describe 
as filamentary crinkles-like structures, roundish flower- and 
ribbon-like features and, at the smallest spatial scales, point- 
like brightenings. Photometric properties and magnetic flux 
density are not constant within most of them, both contrast 
and magnetic flux density being in general higher at their 
edges and lower at their centers. The opposite is observed for 
point-like features. All these stru ctures are highly dynamic 
and c a n evolve one in to another (Rouppe van d er Voort et al. 
12001 IViticchie et ail l2009h . In quiet Sun regions, most 
of the features are point-like and are often clustered in 
groups o r "patches" at the edg es of granules and meso- 
granu les dSanchez Almeida et all [2004 Ide Wijn et alJ I2005L 
2008). The properties of G-band bright features (e.g., dimen- 
sion and contrast) strongly depend on the angular resolution 
of observations so that it is hard to rigorously define the "real 
bright point". For this reason, in the following we will use the 
name "bright feature" to indicate a local excess of brighten- 
ing in the intergranular lanes of G-band filtergrams, indepen- 
dently of the properties of the feature itself. 

Theoretically, brigh t features were investigated through the 
model of flux-tube ( Spruii [T976l) . a small-size concentration 
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of magnetic field surrounded by field-free plasma. Due to 
the high magnetic pressure inside the kG flux tube, the gas 
pressure (gas density) is strongly reduced causing the pho- 
tons to escape from deep photospheric layers that are hotter 
than the surrounding plasma in cold intergranu lar lanes. Nu- 
merical static and dynamic models of isolated dDeinzer et alJ 
1984 t iKnoelker et all [l988t iPizzo et al.1 [19931: ISteiner et al l 
200 U (steinerl 120051) or clustered (ICaccin & Seyerinolll979l: 
Okunev & Kneed 120051: ICriscuoli & Rastl 120091) flux tubes, 
have proved to be successful in reproducing many of the ob- 
served properties, such as the enhancement of contrast toward 
the limb and in some molecular bands (in particular CH and 
CN bands, e.g.. ISchussler et aT] 120031) . the presence of dark 
lanes which often accompany the bright features in off-disk 
center observations, and the suppression of convective mo- 
tions within them. 

From the above figure follows that G-band bright features 
are proxies for kG fields in the solar photosphere. On the other 
hand, the inspection of movies reveals that the brightness of 
the features is not constant with time but intermittently can 
vary such that features can appear dark or bright at different 
instant. This evidence questions the flux tube evacuation as 
the principal physical mechanism responsible for contrast en- 
hancement in small size magnetic concentrations. 

The understanding of the relation between the contrast mea- 
sured in molecular bands and the properties of magnetic con- 
centrations is of relevance for all those studies which em- 
ploy broad-band im aging as proxy of strong magnetic fields 
(de Wijn et al. 2005). In particular G-band imaging of the so- 
lar photosphere, thanks to the high spatial resolution that can 
be achieved in these data, is an independent and complemen- 
tary tool to spectropolarimetric observations. This is of great 
relevance, for instance, for the investigation of the presence 
of kG fields in th e quiet Sun, which i s still a debated topic in 
solar physics (seelde Wr j n et alJl2009l for a review). 

In lViticchie et al.1 (120091) the first r esults from observations 
performed with IBIS dCavallini 2006) in spectropolarimetric 
mode revealed the advantages of combining high spatial res- 
olution spectropolarimetry simultaneous and co-spatial to G- 
band imaging in the study of the temporal evolution of bright 
features in the quiet Sun. In the present work we report further 
results obtained from the same analysis concerning the rela- 
tion between magnetic filling factor and brightness of G-band 
bright features. More in det ail, these comp l ete an d corrob- 
orate the picture outlined in Vitic chie et al.l (120091) only for 
three peculiar cases. The aim of this paper is to go more into 
the details of the capabilities of IBIS and to confirm the find- 
ing of the relation between kG magnetic filling factor and G- 
band bright feature intensity by exploiting the whole analyzed 
dataset. Such a relation will eventually allow us to conclude 
on the fine structuring of G-band bright features. 

In §|2] the dataset is presented, in §[3] the adopted analysis 
methods are described, in § H]we report the results from the 
analysis that are then discussed in §|5j in §|6]the main conclu- 
sions are summarized. 

2. DATA 

A quiet region at disk center was observed on November 2 1 , 
2006 from 16:24 UT to 17: 17 UT with IBIS, at the NSO/Dunn 
Solar Telescope. The field-of-view (FOV) was approxima- 
tively 40" x 40". The acquired dataset consists of 50 scans 
of the two Fe I 630 nm lines performed with a cadence of 
89 seconds. The seeing during the acquisition of the first 36 
scans (which have been analyzed in this work) was excellent 
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FIG. 1. — Two snapshots from the analyzed dataset. Upper panel: G-band 
snapshot; bright feature intensity has been enhanced to show the result of 
the automatic identification procedure. Lower panel: COG magnetogram 
saturated at ±200 G. The images refer to 16:42:29 UT. 

and s table allowing the Adaptive Optics system (Rimmele 
120041) to achieve near diffraction-limited performance. The 
lines were sampled at 45 wavelength points with a spectral 
FWHM of 2 pm and a step of 2.3 pm. In the scanning proce- 
dure the telluric line in between the two Fe I lines was skipped 
while the O2 630.28 nm one was sampled with three points; 
this line allows us to set the absolute wavelength scale of our 
spectropolarimetry measurements. 

To perform spectropolarimetric measurements, the incom- 
ing light to IBIS is modulated by a pair of nematic liquid crys- 
tal variable retarders placed in a collimated beam in front of 
the field stop of the instrument. The light is analyzed by a 
beam splitter in front of the detector, imaging two orthogonal 
states onto the same chip thus allowing for dual-beam spec- 
tropolarimetry. The modulation is in such a way that at each 
wavelength position six modulation states I + S (and its or- 
thogonal states I-S) are detected with the following temporal 
scheme: S = [+V,-V,+Q,-Q,+U,-U]. Spectropolarimetric 
images have a pixel scale of 0.18". Simultaneous to narrow- 
band spectropolarimetric data, broad-band (633.32 ±5 nm) 
counterparts of the same FOV were acquired. The exposure 
time for both narrow-band and broad-band data is 80 ms. 
The pixel scale of broad band data is 0.09". These refer- 
ence images are normally used to register their spectropo- 
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larimetric counterparts. In this particular observation setup, a 
post proc essing Multi-Frame Blind Deconvolution procedure 
(MFBD; Lofda hll2002l) has been applied to these broad-band 
images to obtain a master frame from each scan with much 
reduced seeing degradation and a homogeneous resolution in 
the whole FOV. We are then able to correct for anisoplanatism 
applying a de-stretch process on the single broad-band im- 
ages, using this master frame as reference image. The com- 
puted de-stretch matrixes are successively applied also to the 
spectropolarimetric images. 

Due to the classical mounting of the double etalon sys- 
tem, the spectral bandpass of IBIS experiences a wavelength- 
dependent blueshift within the FOV. For each pixel in the 
FOV, the blueshift amount is computed fitting a template line 
profile on the calibration dataset of flat-fields acquired just 
after the observations, taking into consideration also the dis- 
placements introduced by the de-stretch procedure. The cali- 
bration data includes flats, darks, resolution targets, grids for 
alignment and instrumental polarization calibration measure- 
ments. The IBIS spectropolarimetric I±S images were then 
reduced applying the pipeline developed at NSO and slightly 
modified for this dataset needs, which takes care of dark sub- 
traction, pixelwise gain calibration, co-alignment with broad- 
band images, blueshift correction, de-stretching and also the 
co-registration of all images in each scan and for the whole 
dataset. The frames were then combined and corrected for in- 
strumental polarization to retrieve the Stokes vectors. At this 
stage, also t he telescope induc ed polarization is taken into ac- 
count (as in Jud ge et al.ll20 10). We used the closest DST tele- 
scope calibration, which was measured on February 16, 2007. 

Additionally to broad-band and narrow-band data, G-band 
filtergrams (430.5 ± 0.5 nm) were acquired for a slightly 
smaller FOV; the pixel scale of such filtergrams is 0.037", 
while the exposure time is 15 ms. The fine alignment between 
G-band and broad-band images was performed through grid 
line targets. Tracking inaccuracies were corrected via a cor- 
relation procedure. A post processing MFBD procedure was 
applied also to the G-band images. 

As already stated, the quality of spectropolarimetric data 
has been improved by applying on each spectropolarimetric 
image the same shifts necessary to align and destretch the as- 
sociated broad-band image with respect to the MFBD restored 
broadband image. In this way, the seeing-induced crosstalk is 
reduced and a spatial resolution comparable with that of the 
individual narrow band images is reached. The angular reso- 
lution of the spectropolarimetric dataset has been estimated to 
be 0.4" ~ 2 pixels. 

The average noise level for Stokes V, measured as stan- 
dard deviation of the circular polarization in continuum wave- 
lengths, is oy = 3 • 10~ 3 /c (here Iq is the continuum intensity). 

3. DATA ANALYSIS 

Stokes V measurements are exploited in two different 
ways to conclude on the magnetic properties of G-band 
bright features. Us ing the center-of-gravity method (COG, 
iRees & Semell979l) . we calculate longitudinal magnetic flux 
density maps; the COG method is not affected by saturation in 
the kG regime, for this reason it results suitable for the analy- 
sis of kG bright features. 

Beside the COG approach, the inversion of the Stokes 
/ and V profiles of the two lin es with the SIR code 
dRuiz Cobo & del Toro Iniestalll992l) allows us to determine 
the line-of-sight (LOS) field strength and the magnetic filling 
factor. The polarization profiles emerging from the pixels un- 



der examination are interpreted by means of two atmospheric 
components: one is magnetized while the other is field-free. 
The portion of each pixel occupied by the magnetized compo- 
nent (i.e., the magnetic filling factor f) is a free parameter of 
the inversion. The temperature stratification of each compo- 
nent is modified with two nodes; as initial g uess the Harvard 
Smith sonian Reference Atmosphere model (Gingeri ch et alj 
119711) is used. The LOS velocities and the field strength in 
the magnetized component are both assumed to be constant 
with height. The stray-light contamination is considered to 
be unpolarized and is defined by averaging the Stokes / spec- 
tra in a region of 1 arcsec around each analyzed pixel. Such 
a contribution is added to the Stokes / profile synthesized 
from the model atmosphere weighted by the stray-light factor 
a; the most probable value for the latter factor is a ~ 85%. 
No macroturbulence and microturbulence are considered. Fi- 
nally, the finite spectral resolution of the instrument is taken 
into account using the kno wn spectral PSF of the instrument 
dReardon & Cavallinill2008l) . 

Our analysis is focused on the properties of magnetic bright 
features. For this reason we analyze exclusively very strong 
polarization signals, i.e., Stokes V signals above 4 x ay in 
both Fe I lines. In spite of the severe selection rule, a total 
of about 53000 pixels are considered; these correspond to ap- 
proximately 3% of the whole FOV at any step of the time 
sequence and are found to be in strong correlation with bright 
features in simultaneous and cospatial G-band frames. On the 
other hand, linear polarization signals are found to be always 
below the selecti on threshold so tha t no full Stokes inversion 
is performed. In lBeck et al. ( 2007) the authors were able to 
perform full Stokes inversions in bright feature locations with 
a polarimetric sensitivity three times better than ours for the 
visible lines. 

To automatically identify the bright features in G-band 
filterg rams we adopted the algorithm already used by 
Sanchez Almei da et al.l ((2007) for the identification of bright 
features in 0.2 arcsec angular resolution G-band images filter- 
grams acquired at the DOT. The details of the procedure are 
explained in the paper cited above. 

4. RESULTS 

We investigated the magnetic properties of G-band bright 
features through both COG magnetograms and inversion of 
Stokes V profiles. In the following we present the results de- 
rived from such analyses in two separated sections. The first 
one (§ 14.11) . concerning the analysis of magnetograms, allows 
us a straightforward comparison with the results presented in 
th e ext ended literature on this topic (see §[T|). The second one 
(§ 14.21 ). concerning the inversion analysis, offers new insight 
on the relation between photometric and magnetic properties 
o f G- band bright features sheding new light on the results in 

§EI 

4. 1 . G-band intensity and COG maps 

Fig. Q] (upper panel) shows an example of bright features 
identified through the automatic procedure described in § [3] 
on one of the snapshots of the sequence. The average frac- 
tion of solar photosphere occupied by such features during the 
time sequence is about 0.6%. Their linear dimensions range 
between ~ 0.1 arcsec up to 2-3 arcsec. They appear in a 
great variety of shapes: from the point-like (at the angular 
resolution of the G-band data) to the extended and elongated 
(e.g., the feature in the center of the FOV in the upper panel 
of Fig. ID- 
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FIG. 2._ — Average properties of bright features from both G-band images and COG magnetograms. The plots are derived by calculating: the average flux 
density B mag , the average G-band intensity / (normalized to the average intensity of the time sequence), and the total area A of bright features identified through 
the automatic procedure described in §|5] In the left panel: 7 vs. A scatter plot. In the central panel: 7 vs. B mag scatter plot. In both the plots an exponential fit of 
the data is reported (solid line). The error bars represent the average standard deviation error of the fit; the horizontal dash ed lin e represent the saturation value 
I sat ; the shaded area in the left panel represent the standard deviation error of I ml . For details on the fit procedure refer to jj 14. H and Table[T] In the right panels: 
histogram for / (upper panel) and B mag (lower panel), respectively. 



TABLE 1 

RESULTS FROM THE EXPONENTIAL FITS 



scatter plot a be 

Jvs.A -0.37±0.13 -14 ±7 [arcsec^ 2 ] 1.4±0.1 

I_vs.B_ mag -0.40 ±0.07 -3.0 ±0.9 [10~ 3 G~'] 

/vs./ ' -0.37±0.04 -7.0±0.9 



In Fig. [2] (left panel) we report the scatter plot relating the 
average G-band intensity of bright features / (normalized to 
the average intensity of the time sequence) and their area A. 
The apparent linear relation between I and A is broken for 
A > 0.1 arcsec 2 and / seems to saturate for larger A values. 
Hence, we fitted an exponential function a\ ■ e blX + cy on the 
scatter plot, where a\, b\, and c\ are free parameters of the 
fit. The exponential fit is overplotted on the data while the 
derived parameters are reported in Table [TJ According to the 
exponential fit, we can define a saturation level for the average 
intensity of large bright features approximatively at c\ = T sa , = 
1 .4. In the plot T sa , is represented as a horizontal dashed line, 
while the shaded area represents the error on I sat derived from 
the fit procedure. 

In Fig. |2] (central panel) we report a scatter plot relating 
the average COG magnetic flux density B mag and / of identi- 
fied bright features. In order to compare these quantities, we 
rescaled the COG magnetograms (Fig. [TJ lower panel) at the 
G-band pixel scale, so that we could associate to each identi- 
fied feature its B mag and / computed over the pixels forming 
the region itself. We find that also the relation between B mag 
and / can be described by an exponential law. The fit dis- 
played in Fig. |2] (central panel) was obtained by imposing the 
saturation level to be I sat (i.e., 02 ■ e +Tsat)- The values for ai 
and Z?2 derived from the fit are reported in Table [TJ 

In the rightmost panels of Fig. [2] we report the histograms 
for B mag and I. The most probable value for the average flux 
density in bright features is approximatively 100- 150 G with 
a tail up to ~ 800 G, while the normalized average intensity 
ranges between 0.9 and 1.4 with the most probable value at 
~ 1.1. 

Both the scatter plots in Fig.[2]have been defined by consid- 
ering the data from the whole time sequence. It is important 
to specify that both plots are good representations of the rela- 



tions between the magnetic flux density, the G-band intensity, 
and the area of bright features at any instant of the time se- 
quence. This has been checked by comparing the behavior of 
all 36 scatter plots taken individually along the time sequence 
with the scatter plots of Fig. [2] In Fig. we show one of 
the results of such a check; we report the I vs. A scatter plot 
for each instant of the time sequence in comparison with the 
results of the exponential fit in Fig. [2] (left panel). We have 
to specify that in Fig. [3]/ is obtained by normalizing the in- 
tensity to the average G-band intensity of each frame. The 
agreement between the scatter plots and the exponential be- 
havior derived from the time-integrated scatter plot in Fig. [2] 
is more than satisfactory. Once the consistency has been veri- 
fied, we are allowed to use the time integrated scatter plot, as 
it provides a more statistically sound representation. 

4.2. G-band intensity and inversion results 

Fig. |4] shows the temporal evolution of the magnetic prop- 
erties of the largest bright feature in the FOV, as derived from 
the inversion analysis, together with G-band and COG mag- 
netogram images. The feature can be considered a small net- 
work patch, i.e., a region in which the plasma dynamics con- 
centrates kG fields over a region 2-3 arcsec wide. kG fields 
are revealed over the entire region comprising the bright fea- 
ture (orange contour). High magnetic filling factors are found 
in correspondence with the brightest portions of the feature. 
Namely, in these regions / is mostly higher than 15% and 
can reach values up to ~ 35%. Such values are much higher 
than typical / values for kG fields in the inter network which 
have been measured to be of few percents (e.g., Viticc hie et al.1 
2010). This correlation is persistent along the whole time se- 
quence (cfr., the different snapshots of Fig. |4j>. 

Fig Hshows the pixel-by-pixel correlation between Stokes 
V amplitude, magnetic filling factor /, and G-band intensity 
/ (normalized to the average intensity of the whole time se- 
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FIG. 3. — /vs. A scatter plots along the whole acquired time sequence; 7 is 
normalized to the average G-band intensity of each frame. The time sequence 
is divided in three columns representing three subsequent time intervals; in 
each column the time rolls from up to down. Each scatter plot relates / and 
A for each instant of the time sequence. In all the plots the result of the 
exponential fit shown in the left panel of Fig.f2]is reported (solid line). The 
shaded area here represents the average standard deviation with respect to the 
exponential fit (error bars in Fig. [2). The horizontal dashed line represents 
the saturation value I sat . 



quence) along slices extracted from the analyzed dataset. This 
repres entation is similar to the one adopted by Berg er etaO 
(2004) to represent results from observations performed in 
different bands. Stokes V profiles from both observations and 
the inversion analysis for selected positions along the slices 
are also represented in the figure. These profiles are well fit- 
ted by the SIR code under the adopted inversion hypotheses 
(§ EJ. In the two panels of the figure the variation of the G- 
band intensity is found to be in correlation with both Stokes 



V amplitude and /, while the magnetic field strength does not 
show any correlation with these quantities; in fact, it slightly 
varies around 1.5 kG. The two cases reported in Fig. [5] are 
representative of a point-like bright feature in an intergranu- 
lar lane (upper panel), and an extended bright feature in an 
intergranular lane (lower panel). 

The upper panel of Fig.[5]displays the different capabilities 
of G-band and spectropolarimetric data in resolving the prop- 
erties of bright magnetic features. In fact, the bright feature 
is found to be ~ 0.15 arcsec wide, i.e., four pixels in G-band 
data, while the polarization signal is spread over a region of 
~ 0.7 arcsec, i.e. four pixel in spectropolarimetric data. This 
is a good example of a magnetic feature whose dimension 
is close to the limit of the angular resolution of spectropo- 
larimetric data (§|2). According to the criterion described in 
§|3] we inverted the profiles of three of these four pixels, i.e., 
positions C, D, and E where we measured approximatively 
1 .5 kG, and a filling factor approximately of 10%. 

In the lower panel of Fig. [5] a horizontal slice is extracted 
in correspondence with a strongly magnetized intergranular 
lane. The increase of brightness of the feature from position 
B to D is correlated with the increase of the magnetic filling 
factor, that reaches values above 30%. The Stokes V ampli- 
tude variation measured in positions B, C and D is also ev- 
idently correlated with the G-band intensity variation found 
along the slice. The plots also show that quite high values for 
both magnetic filling factor and magnetic field intensity are 
found in the dark part of the intergranular lane (positions E 
to F). In particular, we notice that at positions B and E the 
estimated filling factor values are similar (about 20%), while 
the observed G-band intensities are quite different. This dif- 
ference has to be ascribed to the spurious contributions to the 
polarimetric signal in pixel E from the adjacent magnetized 
pixels (the two bright features right at the left and below po- 
sition ~ [3,3] arcsec in the G-band subfield). 

The correlation between the magnetic filling factor and the 
G-band intensity shown in Figs. |4] and [5] is found for all the 
bright features observed along the time sequence and selected 
by the automatic procedure. As a further example, in Fig. [6] 
we report the temporal variation of both brightness and mag- 
netic filling factor in a process of merging of faint bright fea- 
tures in a single high brightness feature. To derive the plots 
in such a figure we focused on the properties of the pixels 
around the strongest Stokes V signal in the selected subfield 
(white square on G-band images). This choice was necessary 
since for the first five instants many distinct bright features 
were recognized by the automatic procedure and this does not 
allow one to define a unique temporal evolution. From the 
plot we can recognize a correlation between the increase of 
brightness and magnetic filling factor. 

In Fig. [7] such a correlation is made explicit by the scatter 
plot relating I and the average filling factor / calculated over 
the features selected by the automatic procedure. Also in this 
scatter plot we report the result of an exponential fit performed 
adopting the function 03 ■ e hiX + I sal . The values of 03 and bj 
derived from the fit are reported in Table Q] 

A simple test of the consistency of the results obtained from 
the fits can be done by comparing the values of £>2 and £3 
(second and third lines of Table Q]). In fact, if we consider 
1.7 kG as the typical field strength value derived from the 
inversions, we find that the values reported in Table Q] satisfy 
Z?2 — bi/il.l x 10 3 G). This relation stems from the fact that 
Z?2 and Z?3 are the inverse of a magnetic flux density and a 
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FIG. 4. — Magnetic properties of a large bright feature along the whole acquired time sequence. The time sequence is divided in three columns representing 
three subsequent time intervals; in each column the time rolls from up to down. For each instant of the sequence two frames are reported: a G-band image detail 
of the bright feature under examination (left frame) and a COG magnetogram saturated to ±200 G (right frame). The magnetic properties as obtained from the 
inversion analysis are reported as contours: the kG field region (orange contour), the region with magnetic filling factor > 15% (sky-blue contour), and the region 
with magnetic filling factor > 30% (green contour). The distance between the major ticks of the frames is of 1.4 arcsec. 



magnetic filling factor, respectively. 

5. DISCUSSION 

The automatic identification algorithm applied on G-band 
filtergrams allows us to estimate that ~ 0.6% of the dataset 
FOV is occupied by bright features. This result is in good 
agreem ent with the value reported by Sanchez Almei da et al.l 
(2004) for the quiet photosphere observed at 0.14 arcsec an- 
gular resolution. 

As expected, the SIR inversion retrieves kG fields where 



G-band bright features are found (e.g., Shelvag et al.ll2004l) . 
This co-spatiality, here presented only for three examples of 
features in the FOV, is comm on to all the bright f eatures in- 
dividuated in the dataset (see iViticchie et al.ll2009l for three 
more examples of small features). 

The spectropolarimetric data have a PSF with a FWHM of 
~ 0.4 arcsec, as can be deduced by the shape of the |V| max- 
imum amplitude showed in the upper panel of Fig. [5] There- 
fore, the low filling factor retrieved by the spectropolarimetric 
inversion can be ascribed to the intrinsic size of the magnetic 
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FIG. 5. — Two examples of slice profiles from the results of the inversion analysis and G-band intensity images. These are representative of two interesting 
cases, namely, a ~ 0.3 arcsec wide isolated bright feature observed at 16:42:29 UT (upper panel), and an extended bright feature in a kG magnetized intergranular 
lane observed at 16:39:31 UT (lower panel, see also Figf4j- Left panel: ad xd G-band intensity detail image to identify the position of the slice here represented 
as a white line (left image). Central panel, upper plot: magnetic filling factor / derived from the inversion analysis (solid black line), the horizontal dotted lines 
represent the levels defining the filling factor contours in Fig.f4](i.e., / = 15% and / = 30% levels, respectively); here we plot also the magnetic field strength 
derived from the inversion (blue line) whose variation refers to the 1.5 kG level represented as a dashed line; the shaded area represent a ±150 G field strength 
variation around the 1.5 kG level. Central panel, lower plot: G-band intensity along the slice normalized to the average intensity of the time sequence (black 
line) and \ V\ maximum amplitude along the slice (red line) properly rescaled to be compared with the G-band intensity profile. The histogram-like representation 
adopted in these plots allows us to highlight the different pixel scales of G-band data (0.037 arcsec) and spectropolarimetric data (0.18 arcsec), respectively. 
Right panel: Fe I 6301.5 A Stokes V profiles from both observation and inversion analysis (right panels, symbols and solid line, respectively) picked at different 
positions along the slice (named as A, B, C, D, E, and F); each profile is represented in the x—y range specified in the plot D (wavelengths are represented as 
deviation with respect to the laboratory wavelength of the line). Positions A,B,C, D, E, and F are reported in the G-band detail image as ticks and on the intensity 
profile with dots. 
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FIG. 6. — Temporal variation of both brightness and magnetic filling factor 
in a process of merging of bright features. Left panel, upper plot: temporal 
variation of the magnetic filling factor. Left panel, lower plot: temporal vari- 
ation of bright feature brightness (normalized to the average intensity of the 
time sequence). In both the plots the dotted line represents the trend of the 
evolution obtained by smoothing the temporal variation represented with the 
solid line. In the right panel the time sequence is divided in three columns 
representing three subsequent time intervals; in each column the time rolls 
from up to down. Here, G-band image details of the bright feature under 
examination are reported. The square indicates the region from which the 
quantities in the left panel are derived; this is placed where the strongest 
Stokes V signal is found in the region under examination. At UT 16:54:21 no 
bright features were found by the automatic procedure in the selected square. 



counterparts of the G-band bright features and/or to the PSF 
shape. Even considering the several effects which can con- 
tribute to a PSF with extended tails, namely, the part of the in- 
coming aberrations not corrected by the AO, scattered light in 
the optical path or by the atmosphere, and defects in the image 
calibration pipeline, it seems unreasonable that the PSF effect 
is the dominant one, as we should have a PSF with extremely 
extended tails. Moreover, for the case represented in the upper 
panel of Fig. [5] the filling factors found are compatible with a 
magnetic feature almost filling a single IBIS pixel (0.18 arc- 




0.4 



FIG. 7. — Correlation between average magnetic filling factor / and average 
G-band intensity / (normalized to the average intensity of the time sequence) 
for bright features in the analyzed dataset. In the plot an exponential fit of 
the data is reported (solid line). The error bars represent the average standard 
deviation error of the fit; the horizontal dashed lin e rep r esent the saturation 
value I sa i- For details on the fit procedure refer to § 14. ll § 14.21 and Table [TJ 



sec) spread by a 0.4 arcsec FWHM wide PSF. Therefore, we 
interpret the ~ 15-30% filling factors found for the extended 
bright feature (Fig. |4]and lower part of Fig. |5j as a signature of 
unresolved magnetic features. A deeper analysis of the PSF 
convolution effects on spectropolarimetric inversions will be 
the topic of a forthcoming paper. 
From the unsigned flux density histogram, we can infer that 
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the most probable unsigned flux density of the bright fea- 
tures is approximatively 100- 150 G (Fig. [2] ri ght panels). 
This v alue is about twice the v alue obtained by iBeck et al.l 
(2007), but is in agreement with iBerger & Titld (120011) . The 
disagreement among the results can stem from the differ- 
ent spatial resolution of the magnetograms: the spectropo- 
larimetric data employed in IBeck et al] (120071) had a spatial 
re solution of about 1 ar c sec, while the magnetograms used 
in IBerger & Titld d2001l) had a 0.3 arcsec spatial resolution, 
much closer to that of our dataset. 

When plotting the G-band intensity vs. the magnetic flux 
density (Fig. central panel), we find that the G-band in- 
tensity seems to saturate for values of B mag > 500 G. We 
can compare this plot with Fig. 12 of [Berger & Titld (1200 ll) 
and the G-band contrast of identified magnetic brightening 
as a function of th e peak magnetic flux density in Fig. 9 of 
IBerger et al.l (120071) . with the caveat that the authors consid- 
ered peak quantities instead of average quantities. The scat- 
ter plot relating the G-ba nd contrast and the total integrated 
polarization reported bv lBeck et al.l (120071) also shows an ev- 
ident saturation of the contrast around 0.53 for the reported 
values of the total polarization. In those plots, even if the 
contrast values are in agreement with the ones we measured, 
the authors did not find any relation between the two quanti- 
ties. We can note, howev er, that our statistics is much richer 
than in lBeck et al.l (120071) . Such a difference in the statistics 
stems from two main facts. The first is that we analyzed 35 
sets of co-spatial an d co-temporal G-b and and spectropolari- 
metric data, while in Beck et al] (120071) 8 sets were available. 
Moreover, different pr ocedures for the se lection of bright fea- 
tures were adopted; in IBeck et aT] (120071) the authors specify 
that their procedure usually leads to the selection of extended 
patches that could be divided in smaller ones, this caused an 
important reduction of the total number of selected bright fea- 
tures. 

In particular, the latter two works lack statistics exactly for 
low magnetic flux densities and try to fit their scatter plots 
with linear functions, while the former, on the contrary, only 
shows magnetic flux densities below 400 G. 

Commenting the plot in Fig. [2] (left panel), the in- 
crease of G-band intensity with the increase of the size 
of bright features is in general agreement with previ- 
ous results from high spatial resolution observations (e.g. , 
Berger et al.lll995HBovelet & Wieldl2003t IWiehr et al.l l2004t 
Hirzberger & Wiehri 120051: IBerger et all 120071) . The G-band 
contrast saturation values reported are usually between 1.2 
and 1 .4 for disk-center observations, nevertheless, due to the 
large scatter of data, dif ferent interpretatio ns of the results 
have b een suggested. In Wie hr et aiT(l2004l ) and IBerger et alj 
( 1995) for instance, the increase was within the dispersion of 
the results and the authors concluded that no clear trend of in- 
tensity as a function of size was observe d, the contrast being 
essentially constant. Ber ger et"afl ((2007) sho wed instead that 
their r esults can be fitted by a straight line. Bove let & Wiehrl 
(2003) reported a clear increase of contrast with size and fit- 
ted the results with a polynomial fun ction of third degree. The 
results presented in lUtz et ai] d2009al) are an exception in this 
context, since the authors reported a decrease of contrast for 
bright features larger than HINODE resolution. The authors 
interpreted this trend as an effect of the segmentation algo- 
rithm employed to detect the analyzed features. Close in- 
spection ofjdlthese data (with the exception of the ones in 
Utz et al. 2009a) reveals an increase of contrast when the size 
of the feature is comparable with the dataset resolution and 



a saturation at the largest sizes. This is in agreement with 
our finding and we can notice that even the threshold spa- 
tial scale value of contrast saturation that we obtain from our 
analyses is compatible with previously presented results. In 
fact, assuming radial symmetry for the selected features, we 
find a threshold size value of approxima tely 2x^/1 /(irb\) ~ 
0.3 arcsec ~ 210 k m, in agreement with IBerger et all d 1995b . 
IWiehr et al] d2004l) . and lBerger et al.l d2007l) . when taking into 
account the different spatial resolution of the data. 

This contrast/size relation is in apparent disagreement 
with classical flux tube models, according to which con- 
trast decreases with the size of the magnetic fe ature (e.g., 
Spruit! 1 9761; iDeinzer et ai]ll984l;lPizzo et al.lll993l) . Recently, 
Criscuoli & Rastl d2009i) showed that this effect can be inter- 
preted as a signature of unresolved magnetic features. The au- 
thors, using a two-dimension numerical model of isolated and 
clustered magnetic flux tubes, simulated the emergent inten- 
sity for different cluster dimensions at different angular reso- 
lutions. They found an exponential like relationship between 
size and continuum contrast, which saturated for cluster di- 
ameters larger than twice the spatial resolution of the obser- 
vations. This picture is also corroborated by the fact that the 
finite spatial resoluti on of observations most ly affects small 
size features (e.g.. ICriscuoli & Ermolli 2008, and references 
therein): the smaller the feature the larger is the reduction 
of its contrast with the decrease of spatial resolution. Larger 
features are instead less affected, thus explaining the general 
agreement of contrast values presented in the literature at the 
saturation. 

The work of IWiehr efafl d2004l) . which showed that both 
the 588 nm continuum contrast and the G-band contrast of 
bright features saturate at the sa me spatial scale, allows us to 
directly compare our result with ICriscuoli & Rast (20Q9|) and 
the two results seem to match quite congruously. In fact, we 
found the saturation for A > 1/bi ~ 0.07 arcsec 2 (Table[T|), for 
which we can calculate an equivalent diameter of d ~ 0.3 arc- 
sec, i.e. about twice the spatial resolution of our G-band 
dataset. 

The plot reported in Fig. [7] showing the relation between 
magnetic filling factor and G-band intensity found in our data, 
also supports the idea of the presence of unresolved mag- 
netic features in our observations. We remember that, only 
the Stokes profiles cospatial with pixels which show G-band 
brightness enhancement have been considered to define the 
plot and that the angular resolution of the spectropolarimetric 
data has been estimated to be ~ 0.4 arcsec. From these inver- 
sions, it follows that filling factor values are mostly < 10% 
and only a small fraction of the bright features present filling 
factors above 20%. Beside this, the G-band filtergrams have 
a spatial resolution of ~ 0.1 arcsec, therefore in these images 
we can resolve bright features that are as small as 25% of 
the spectropolarimetric data pixel. The low values of the fill- 
ing factor suggest that the smallest scale bright features are 
smaller than 0. 1 arcsec. 

The variation of the G-band intensity for point-like features 
(A < 0.1 arcsec 2 ) supports this idea: fully resolved bright fea- 
tures should not change their intensity with their dimension. 
We recall that, in G-band bright features, the brightness en- 
hancement occurs because of the local evacuation induced by 
the high magnetic pressure. We can assume that increasing 
the magnetic filling factor means enlarging the evacuated re- 
gion in our pixel and therefore interpret the plots of Figs. 01 
|6l|5]and|7]with the consequent enhancement of the number of 
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photons that emerge from hot deep layers in the photosphere 
in the pixels forming the bright feature under examination. 
These arguments, combined with the G-band intensity satu- 
ration found for A > 0.1 arcsec 2 and / > 20%, suggest that 
the larger G-band bright features can be thought as clusters 
of kG "elementary bright points" having ty pical dimension 
< 0.1 arcsec. This is in full agreement with Criscuoli & Rast 
(200§. 

From our inversion analysis we conclude on the role of the 
magnetic filling factor of kG fields in the formation of G-band 
bright features (Fig. |7): the larger the magnetic filling factor 
the brighte r the feat u re. T o explain our findings, we can re- 
fer to Ide Wijn et al.1 (120051) . The authors commented about 
the observation of G-band bright features in quiet Sun (inter- 
network) regions: "we conclude that, even though the inter- 
network bright points may momentarily become invisible, the 
magnetic field element remains and may become bright again 
at some later t i me. T his agrees well with the conclusion of 
\Berger & Title\ < \2001\) that magnetism is a necessary but not 
a sufficient condition for the formation of a network bright 
point". The results here presented in § [4] complement their 
argumentation. In terms of the spatial extension of kG flux 
tube clusters the sentence can be reworded as: the larger the 
cluster the brighter it will appear, but, once the cluster dimen- 
sions are comparable with the spatial resolution of the obser- 
vations, the brightness saturates. From this follows that small 
G-band bright features can be momentarily invisible simply 
because the local spatial organization of kG fields is temporar- 
ily changed, for example, through the action of the photo- 
spheric dynamics or because the spatial resolution decreases 
due to seeing. This affects the stron gest the brightness o f 
point-like features (as already shown in Vitic chie et al.l200"9l) . 
but also can induce the internal bright ness variation of large 
feature s (Figs. [4] and |5}. In the recent Sanch ez Almeida et al.l 
(12010 ) the measure of bright feature abundance in the quiet 
Sun is based on the concept that small scale faint features can 
be extremely variable in time. Applying an identification pro- 
cedure that makes use of 0. 1 arcsec angular resolution G-band 
filtergrams time sequences, the authors were able to rise the 
estimate of G-band bright feature abundance of a factor three, 
i.e., up to ~ 1 BPs every 2 arcsec 2 . 

It is worth also to focus on recent alternative interpretations 
of the magnetic structure of photospheric bright features. In 
the last decade many works based on modern MHD simula- 
tions performed with the MURaM code have been dedicated 
to the study of G-band bright features ( [Schussler et alj|2003t 
IShelvag et al.ll200l: iVogler et al.ll200l iShelvag et alJl2007l) . 
These were able to reproduce the photometric properties of 
G-band observations and to confirm the association between 
local brightness enhancement and small scale kG concentra- 
tions. In all these works a common figure for the magnetic 
structure of bright features emerges: at 20 km spatial reso- 
lution these are "elongated sheet-like" features. No evidence 
of internal structuring, e.g., through many point -like features 
is foun d. In agreement with this interpretation, Berger et al. 
(2004) reported that no evidence of internal structuring was 
found for bright features analyzed in both 0.1" resolution G- 
band images and 0.18" resolution magnetograms obtained 
at SST In contrast with our interpretation, the authors put 
forward that an internal variation of brightness of extended 
bright features can be explained via an internal variation of 
field strength instead of magnetic filling factor. In the re- 
cent work of Narayan & Scharmer] (120101) the authors ana- 



lyzed full Stokes measurements performed through CRISP at 
Solar Swedish Tower (SST ) with 0.15 arcsec a ngular reso- 
lution. On the same line of Berge retafl (120041) . the authors 
analyzed Stokes V profiles adopting a Milne-Eddington code 
and by imposing the magnetic filling factor to be equal to 
unity, i.e., considering magnetic flux concentrations to be re- 
solved in their observations. In their work the authors prop- 
erly discussed this choice specifying that their measure of the 
magnetic field should be consider ed either as a me asure of 
the magnetic flux density (as in Berg er et al.l 120041 ) or as a 
measure of a roughly defined average of the magnetic field 
strength. On the other hand, our approach considers the ef- 
fect of the limited spatial resolution in the inversion analysis 
of 0.4 arcsec spatial resolution data; this allows us to separate 
the contribution of the field strength and the magnetic filling 
factor. 

The work s of iNaravan & Scharmer (120101) . based on 
CRISP data, iDanilovic et al.l (12010a!) and other observations 
performed with SUNRISE/IMaX, represent the first effort in 
trying to understand the photospheric magnetism at 0.15" an- 
gular resolution. Such observations will shed new light in the 
study of the structuring of the photospheric magnetic field. 

As a conclusive part of our discussion we want to point out 
two recent works which support a fine structuring of photo- 
spheric bright features. In the first one, iGoode et all 12010), 
the observation of point-like bright features in 0.12" reso- 
lution observations at 705.7 nm is reported, while no sheet- 
like features were obse rved. The second one is the work of 
IDanilovic et al.1 (1201 Obi) in which processes of magnetic inten- 
sification are simulated via the MURaM code. Such processes 
form stable point-like (i.e., ~ 0.1" wide) kG features which 
are able to produce a continuum brightness enhancement at 
630 nm. These are different from the strong sheet-like kG 
concentratio ns produced by advert ing field lines in downflow 
regions (e.g., Schuss ler et al.l2003l) . From our point of view, it 
could be extremely interesting an in-depth study of the spatial 
and temporal evolution of such kG elements to check whether 
they can produce, by grouping, extended ~ 1" wide bright 
features. 

6. CONCLUSIONS 

We have presented new results from the analysis of observa- 
tions with IBIS in spectropolarimetric mode. These are, at the 
same tim e, complementary and corroborative of the results re- 
ported in Viticc hie et alj (2009). Namely, in the present work 
we confirm the correlation between bright feature brightness 
and kG magnetic filling factor found in the analysis of three 
p articular cases of tem poral evolution of small-scale features 
in lVitic chie et al] d2009l) : the higher the local kG filling factor 
the brighter the G-band feature. This conclusion, first evinced 
by the detailed analysis of an extended bright feature (i.e., a 
small network patch), has been validated by a statistical study 
of the properties of all the bright features observed in the FOV 
at any instant of the time sequence. 

Considering the angular resolution of our spectropolarimet- 
ric (~ 0.4 arcsec) and G-band (~ 0.1 arcsec) data, we put 
forward the upper limit of ~ 0. 1 arcsec for the typical dimen- 
sion over which the "elementary G-band bright features" (i.e., 
bright points) are formed. Adopting this figure, larger G-band 
bright features could be thought to be cluste rs of bright points 
(in agreement with ICriscuoli & Ras t 2009). Such substruc- 
turing of bright features can be just inferred from 0.4 arcsec 
angular resolution IBIS spectropolarimetric data. New higher 
angular resolution observations will offer to the solar commu- 
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nity the opportunity to conclude on the magnetic structure of 
G-band bright features. 
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